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Contemporary models of hadronic interactions are reviewed. Basic phenomenological 
approaches are compared, with an emphasizes on the predicted air shower characteristics. 
Special attention is payed to the remaining discrepancies between present hadronic MC 
generators and cosmic ray data. Finally, future prospects concerning model improvements 
are discussed, in particular, regarding the possibilities to discriminate between different 
' models on the basis of accelerator or cosmic ray measurements. 

<N ' 

£>. ■ 1 Introduction 

o 

£f~) . Nowadays Monte Carlo (MC) models of hadronic interactions are extensively 

s ! ' used in collider and cosmic ray (CR) physics, being applied for designing new 

experiments, analyzing and interpreting data, or for testing new theoretical ideas, 
when the latter are incorporated in the corresponding MC generators. Here the 
last option is more characteristic for accelerator physics. High energy cosmic rays 
are typically studied by means of extensive air shower (EAS) techniques: instead of 
detecting primary CR particles directly one measures different characteristics of the 
nuclear-electro-magnetic cascade, so-called air shower, induced by their interactions 
in the atmosphere. Because of a large number of hadronic collisions during such 
a cascade it is very difficult to disentangle partial contributions of interactions of 
particular energies and to test individual features of MC models. In CR physics 
hadronic generators are an important technical tool, which extrapolates current 
theoretical and experimental knowledge from the energies of present accelerators up 
to the highest CR energies and helps to establish a relation between the properties 
of primary cosmic ray particles and the observed EAS characteristics. 

Are there significant differences between collider and cosmic ray models? Is it 
possible to apply popular accelerator codes for CR data analysis? As we shall discuss 
below, all contemporary hadronic MC generators are based on similar physics ideas. 
Both at colliders and in cosmic rays one employs various approaches to describe 
hadronic and nuclear interactions; in neither case there exist a preferred approach, 
inherent for this or that particular field. However, generally accelerator codes are 
not suitable for an immediate application in CR field, being unable to treat some 
necessary reactions, like interactions of pions and kaons with nuclei or minimum bias 
nucleus-nucleus collisions, to perform predictive calculations of total inelastic and 
diffractive cross sections, or to be extrapolated into the ultra-high energy domain. 
It is worth stressing again that cosmic rays experiments can test model validity 
only indirectly. Analyzing CR data one has to rely on model predictions, without 
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real possibilities to re-tune model parameters or to refine corresponding algorithms. 

The goal of this work is to review basic approaches employed in contemporary 
CR models, to analyze remaining problems and to outline promising directions for 
future model improvements and for model discrimination on the basis of accelerator 
and cosmic ray measurements. 



2 High energy interactions: "pure QCD" approach 

Nowadays perturbative quantum chromo-dynamics (pQCD) is the only strict 
theory which provides the basis for both qualitative and quantitative description of 
hadronic production in high energy reactions. However, precise pQCD calculations 
are only possible for a restricted class of processes, namely those which are char- 
acterized by a large momentum transfer Q 2 , such that running coupling a s (Q 2 ) 
becomes small and corresponding contributions can be re-summed to a reasonable 
accuracy pQ. Thus, one can not apply it to general hadron-hadron processes, which 
are expected to be dominated by the "soft" (low Q 2 ) physics. On the other hand, 
in the very high energy limit one expects that the role of high virtuality (high 
p 2 ) processes greatly increases, as their contributions are enhanced both by large 
parton multiplicity and by large logarithmic ratios of transverse and longitudinal 
momenta for successive parton emissions ^El- Hence, one may develop MC gen- 
erators, which rely on the pQCD formalism in the description of "hard" (high p t ) 
processes while treating the "soft" ones in a simplified fashion, assuming that the 
latter play smaller and smaller role when the energy increases. 

Quite often such generators are referred to as "pure QCD models" . How much 
is defined there by the perturbative QCD? Typically the main pQCD input is the 
inclusive cross section for production of parton jet pairs with transverse momenta 
larger than some cutoff Qq, given to leading logarithmic accuracy as 

Jot/ .2. _ \- [dp 2 f dx+dxj da 2 f 2 (x+xj S ,p 2 ) 



x r a (4,p 2 ) fi(xj,p 2 ) e(p?-Oo) 



where s is the cm. energy squared for the interaction, da 2 ~* 2 /dp 2 is the lowest order 
parton-parton scattering cross section, and j % a (x,g 2 ) is the parton i ((anti-)quark 
or gluon) momentum distribution function (PDF) in hadron a, when probed at the 
virtuality scale q 2 . Here partons i and j in Eq. JIJ, participating in the high p t jet 
production, originate from partonic cascades in hadrons a and d, being emitted by 
parent partons of smaller p t and larger energy, those in turn being emitted by their 
own parents and so on. The whole process can be described as a parton ladder, 
which starts in hadrons a and d with partons of p\ ~ Qq. 

Still, to develop a MC model one has to make a number of additional assump- 
tions, which are the essence of the mini-jet approach j3J 2] employed in SIBYLL 
pi El and in a number of popular collider models, like PYTHIA [7j or HIJING [Sj. 
First, as (s, Qq) increases with energy faster than the observed total cross sec- 
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Fig. 1. Hadron-hadron scattering as multiple exchange of parton ladders. 



tion, it appears mandatory to consider multiple jet production processes, in other 
words, to describe the interaction as multiple exchange of parton ladders, see Fig.^ 
Furthermore, one suggests that parton distribution in hadron a is governed by 
its electro-magnetic form factor pl/ m and, as a consequence, that at given impact 
parameter b average number of jet production processes is defined by the product of 
cr^J and the hadron-hadron overlap function A a d(b) = J d 2 b' p e J m (b') p e J m (\b — b'\)- 

Kd{s,b,Ql)) = j*{s,Ql) A ad (b) (2) 

Finally, assuming Poisson distribution for the number of jet processes at given b, 
with the average value (« J a c j(s, b, Qf,)) = 2x|^ rd (s, b), and adding a similarly defined 
contribution of "soft" processes X a °d \ s i = ^c soit (s) A a d(b), with a parameterized 
bare "soft" cross section cr soft (s), one can express total inelastic cross section via 
the eikonal X ad = Xad* + X a T as [21 IE] 

= J d 2 b [l - e- 2 *- W>] (3) 

Using this scheme one can generate inelastic events, starting from sampling 
the impact parameter b for the interaction, defining the number of hard processes, 
and generating parton cascades which preceed the hardest parton-parton scattering 
or which follow after it |S| ■ Then, one employs phenomenological procedures to 
convert final partons into hadrons, as well as to treat the soft part of the interaction. 
The most popular one is the string picture 9 : when final partons move apart, being 
connected to each other by color field, a color string is stretched between them. As 
long as the distance between partons increases so does the string tension, which 
finally gives rise to the break up of the string and to a multiple creation of quark- 
antiquark and diquark-antidiquark pairs, those forming secondary hadrons. 

The disadvantage of the approach is a somewhat artificial separation of hard 
and soft processes and a rather arbitrary treatment for the latter. For the former 
one considers parton cascades as starting at given cutoff scale Qcm thus neglecting 
the contribution of low pt (pt < Qo) partons. 

3 High energy interactions: Reggeon approach 

In Gribov's Reggeon approach JU] hadronic collisions are described as multiple 
scattering processes, where each elementary re-scattering corresponds to a micro- 
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scopic parton cascade, just like in Fig. ^ Such an elementary process is treated 
there phenomenologically, as an exchange of a quasi-particle, the Pomeron, so that 
the scattering amplitude is the sum of contributions of graphs of Fig. O 




Fig. 2. Typical contribution to hadron-hadron scattering amplitude; elementary re- 
scatterings are treated as exchanges of composite objects - Pomerons. 

The Pomeron amplitude is postulated in the form 

flifab) = l7 a7^ Qp (°^ 1 /Aa rf (s)exp[-6 2 /(4A ad ( S ))] EE* x P d ( s ,6) (4) 
X ad {s) = Rl + R 2 d + a' P {Q)\ns, (5) 

which is characterized by a power-like energy rise and by a logarithmically in- 
creasing slope X a d{s)- Here ap(0), a'p(O) are the Pomeron trajectory parameters 
(intercept and slope) and j a , R\ are the parameters for Pomeron- hadron a coupling. 

Having specified the structure of Pomeron-hadron vertices, one can obtain the 
total elastic amplitude for hadron-hadron scattering, hence, via the optical the- 
orem, also total cross section. Applying the Abramovskii-Gribov-Kancheli (AGK) 
cutting rules ^2]> one can separate different contributions to the total cross section, 
corresponding to particular final channels of the reaction. For example, assuming 
cikonal vertices and neglecting diffraction dissociation and inelastic screening, one 
obtains Poisson distribution for the number of elementary production processes at 
given b and arrives to Eq. for inelastic cross section, with Xad{s,b) = Xa*d( s >b)- 

How to match this scheme with the pQCD treatment? In general, parton cas- 
cades contain both "hard" {p\ > Qq) and "soft" (j>\ < Qq) parts, Qq being a chosen 
cutoff for pQCD being applicable. In the "semi-hard Pomeron" scheme, employed in 
QGSJET J2j and NEXUS ^3] models, one applies phenomenological soft Pomeron 
treatment for the soft part of the parton cascade and describes general semi-hard 
processes as exchanges of a "semi-hard Pomeron" , the latter being represented by a 
piece of QCD ladder sandwiched between two soft Pomerons [1311151 ITS] , see Fig- 01 
A general Pomeron appears to be the sum of the semi-hard and soft ones, the latter 
corresponding to a cascade of low pt (pt < Qo) partons. 

The scheme coincides with the mini-jet approach in the treatment of the hard 
sector, represented by the ladder parts of the semi-hard Pomerons, see Fig. [21 An 
important difference comes from considering multiple soft interactions, described 
as soft Pomeron exchanges, and from accounting for an additional source of particle 
production, which comes from the low p t part of the parton cascades ("soft pre- 
evolution") in semi-hard processes, shown as soft Pomeron "blobs" in Fig. [31 There 
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Fig. 3. A general "Pomeron" (l.h.s.) is the sum of the soft and the semi-hard ones - 
correspondingly the 1st and the 2nd graph in the r.h.s. 

is no artificial separation between soft and semi-hard processes, both being partial 
contributions to a general parton cascade, depending whether the latter develops in 
the low p t region or extends to higher virtualities, pf > Qq. The same soft Pomeron 
asymptotics is used to describe both soft processes and soft parts of semi-hard ones, 
defining in the latter case corresponding PDFs at scale Qo |14Mlb) . 

In general, there is no sharp border between the models of mini-jet or semi- 
hard Pomeron type. For example, modern version of SIBYLL jS] takes into account 
multiple soft interactions, described as Pomeron exchanges; DPMJET J7| considers 
an additional soft Pomeron exchange for each semi-hard process, which is not very 
different from the above-described "soft pre-evolution" treatment. 

4 Non-linear effects in high energy collisions 

Treating interactions at very high energies and small impact parameters one 
deals with the regime of high parton densities. There, individual parton cascades 
are no longer independent of each other; their overlap and mutual influence give 
rise to significant modifications of the picture described so far. At microscopic level 
such non-linear effects are described as merging of parton ladders [2] , see Fig. 0] 




Fig. 4. Examples of diagrams which give rise to non-linear parton effects. 

The corresponding pQCD treatment is based on the assumption that par- 
ton densities in the low virtuality region are saturated and any further parton 
branchings are compensated by the fusion of individual parton cascades In 
mini-jet models one typically employs a phenomcnological approach, introducing 
some parameterized energy- dependence of the p t -cutoff for mini-jet production, 
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Qo = Qo(s), the latter being an effective saturation scale [Hj- Still, completely ne- 
glecting the contribution of the saturated region (partons of pt < Qo(s)) to particle 
production may be a too strong assumption. On the other hand, applying the same 
Qo(s) cutoff to all collisions at a given energy appears to be a rather crude ap- 
proximation; one loses correlations with actual parton densities, which depend on 
the parton Feynman x, on the impact parameters, and on the projectile and target 
mass numbers in case of nuclear collisions. 

A dynamical description of non-linear interaction effects has been proposed in 
the QGSJET-II model |18| . based on the assumption that saturation effects can be 
neglected for parton virtualities bigger than some fixed, energy-independent cut- 
off Qo. In this scheme the multi-ladder graphs of Fig. Q] are replaced by enhanced 
diagrams, which correspond to Pomeron-Pomeron interactions |191 120| . Suggesting 
that multi-Pomeron vertices are dominated by parton processes at comparatively 
low virtualities \q 2 \ < Qq, one can re-sum all significant contributions of that kind 
and develop a self-consistent MC generation procedure for hadronic and nuclear 
collisions IB| . The main parameter of the scheme, the triple- Pomeron coupling, 
has been inferred from HERA data on hard diffraction in deep inelastic scattering 
reactions. By construction, non-linear screening corrections appear to be correlated 
with corresponding parton densities and become larger at higher energies, smaller 
impact parameters, or for collisions of heavier nuclei, resulting finally in the satu- 
ration of PDFs at the scale Qq and in a considerable reduction of "soft" particle 
production. In particular, the predictions of the QGSJET-II model stay in agree- 
ment with the data of the RHIC collider on the multiplicity of secondary hadrons 
produced in central nucleus- nucleus interactions |21j . 

5 Model sensitivity of air shower calculations 

Applying different models to air shower calculations one typically obtains some 
spread of the predicted EAS characteristics. To get insight into such differences, one 
compares different models with respect to the interaction characteristics which are 
most relevant for air shower development, i.e. regarding the predicted energy rise 
of inelastic cross sections, the relative energy loss of "leading" secondary particles 
(the inelasticity), or the multiplicity of produced secondaries. Such comparisons 
help to understand discrepancies between calculations and experimental data and 
to outline directions for future model improvements. 

On the other hand, not every problem of EAS data interpretation may be at- 
tributed to potential model deficiencies. Nowadays there exist a bulk of accelerator 
data which seriously constrain predictions of hadronic models. Moreover, due to 
a large number of interactions in a hadronic cascade, observed air shower char- 
acteristics appear to depend rather on the general interaction trend over a wide 
energy range than on a sudden change of the interaction mechanism at a particular 
energy. In particular, comparatively low energy interactions, being well studied at 
accelerators, contribute with a large weight |22j- Thus, it would be naive to expect 
that some new valid model will change basic EAS predictions by a large amount. 
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To quantify these statements one can perform simple tests with presently avail- 
able models. For example, using the QGS JET model with the corresponding proton- 
air cross section being artificially rescaled by ±10%, while keeping pion-air cross 
section unchanged, one obtains, due to the shift of the primary particle interaction 
point, about 10 g/cm 2 variation of the predicted shower maximum depth -A max ; 
the corresponding change of the charged particle number N e at ground ranges from 
6% at 10 16 eV to 2% at 10 19 10 20 eV, as shown in Fig. For charged particle 




Fig. 5. Left: shower maximum depth for proton-induced EAS as calculated with the 
default QGSJET model (full line) and with ±10% modifications of eri^air (dashed lines). 
Right: variations of the calculated electron number at sea level for proton-induced vertical 
EAS for ±10% modifications of (relative to the default QGSJET). 

density at 1000 m distance from the shower core p c h(1000) such a sensitivity is 
even weaker, falling below 1% above 10 18 eV, as the reduction (enhancement) of 
N e at ground is compensated by a flatter (steeper) lateral distribution of particles 
in case of larger (smaller) cross section. Thus, present experimental techniques, 
which reconstruct the primary energy on the basis of p c h(1000), should not exhibit 
any significant model dependence. A similar conclusion holds for the fluorescence 
light-based energy reconstruction procedures |28| . 

Still, any studies of the cosmic ray composition, either based on A max mea- 
surements or making use of electron-muon number correlations, exhibit significant 
model sensitivities. As for the shower muon number at ground level, it is very 
robust with respect to the discussed cross section variation, as shown in Fig.HJleft), 
but depends on the multiplicity of hadronic interactions. One can consider two al- 
ternative modifications of model predictions - either increasing the multiplicity of 
just the primary particle interaction by a factor of two or enhancing the multiplicity 
of all secondary pion-air and kaon-air interactions by 10% - see Fig. fright). As is 
easy to see from the Figure, both changes make a similar impact on the predicted 
EAS muon number at highest energies. Evidently, to obtain a factor of two increase 
of Nfj, one needs a comparable enhancement of multiplicities of both primary and all 
secondary interactions, which would strongly contradict available accelerator data. 

What are characteristic differences between presently available models? QGSJET 
and SYBILL predictions for o^L differ from each other by 10-=- 15% at 10 19 eV; the 
disagreement for the multiplicity of charged particles N£t ail also increases with en- 



Czech. J. Phys. 51 (2001) 



A 



7 



S S Ostapchenko 



s. 1.1 
5* 



o M -10 



n in=1 +10% 



N .+100% 



Fig. 6. Sensitivity of the muon number (E^ > 1 GeV) at sea level for proton-induced 
vertical EAS to the cross section and multiplicity variations (as discussed in the text). 




Fig. 7. Inelastic hadron-nitrogen cross section (left) and multiplicity of charged particles 
(right) for p — 14 A and tt — 14 A collisions as calculated with the QGSJET-II, QGSJET, 
and SIBYLL models - full, dashed, and dotted curves correspondingly. 



ergy, reaching a factor of two at ultra-high energies, as shown in Fig.0 Main sources 
of these differences are the treatment of semi-hard processes in the models, which 
depends on the chosen PDFs, and the implementation of non-linear effects. The 
latter are not taken into account in QGSJET, which employs flat ("pre-HERA") 
parton distributions. SIBYLL uses realistically steep PDFs, which leads to a faster 
energy increase of <7j£!^ r . On the other hand, the faster multiplicity increase in 
QGSJET is due to the semi-hard Pomeron scheme employed, where additional par- 
ticle production comes from the "soft pre-evolution" , i.e. from hadronization of 
partons of small pt (pt < Qo)- It is noteworthy that cross sections in QGSJET- 
II are rather similar to the ones of QGSJET - see Fig. [Tfleft), despite the fact 
that the model uses realistic PDFs compatible with HERA data. This is because 
the effect of new PDFs is compensated by non-linear screening corrections. Such a 
compensation does not hold in SIBYLL, where such corrections are introduced via a 
Qo(s) cutoff for hard processes and are neglected for the "soft" component. On the 
other hand, non-linear effects reduce particle production in QGSJET-II compared 
to QGSJET, bringing it closer to SIBYLL - see Fig. fright). For EAS characteris- 
tics one obtains a rather large difference between QGSJET and SIBYLL predictions 
at highest energies, which reaches in particular 30% for the muon number. In turn, 
for QGSJET-II this difference is reduced to 10% [H]. 
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6 Some remaining problems 

Among the most remarkable recent results are KASCADE studies of electron- 
muon number correlations |24) . which, apart from providing information on the 
spectra and mass composition of cosmic rays in the knee region, allowed one to 
quantify the discrepancies between present MC model predictions and CR data. 
The analysis showed that both QGSJET and SIBYLL have problems with the 
data interpretation. The corresponding discrepancies stay at the level of 10% for 
the predicted N e and and arise for the two models considered at different 
energy intervals. To reach a better agreement with the data it is desirable to have a 
model which predicts slightly deeper Jf max compared to QGSJET, or, alternatively, 
reduces N^, moving towards SIBYLL predictions. In fact, both holds in QGSJET- 
II, where non-linear effects sizably reduce particle production compared to original 
QGSJET. Therefore, one may expect of it a somewhat better performance here. 

Apart from that, present models seem to be unable to describe high multiplicity 
muon bundles observed by LEP collaborations One could think of improving 
the situation via a significant enhancement of multiplicity of hadronic interactions. 
On the other hand, as the highest multiplicity bundles mainly originate from events 
with the shower core inside the detector, one may expect a better performance from 
models which predict a deeper -Y max , hence, a steeper muon lateral distribution 25 . 

Finally, there exist a long-standing contradiction between ground-based and 
fluorescence-based EAS detection techniques, concerning both obtained CR spectra 
and mass composition While the composition issue is still difficult to address, 
the energy reconstruction methods, as discussed above, should not exhibit signif- 
icant model dependences. Hence, further studies of systematics of corresponding 
experiments seem to be desirable. 

7 Future prospects 

During last decade one observed a convergence of EAS predictions of contempo- 
rary cosmic ray interaction models. The remaining differences are still significant, 
resulting in particular in sizable uncertainties of CR composition studies. On the 
other hand, numerous new data, obtained with cosmic rays or at colliders, open the 
way for model discrimination and for their further development. Among the most 
important model uncertainties are those connected to the energy extrapolation of 
inelastic cross sections and of the inelasticity of hadron-air interactions. The new 
measurement of Oj^air by tne HIRES collaboration ,2Z] allows us to reduce the for- 
mer and points towards smaller cross sections than those of present models. Future 
LHC measurements of will give us possibilities to calibrate reliably correspond- 
ing model procedures. On the other hand, RHIC collider provided valuable data 
on baryon "stopping" in central nucleus- nucleus collisions j2Hj- When repeated for 
different centrality selections and for various combinations of projectile and target 
nuclei, those will allow us to fix reliably model predictions for the inelasticity of 
hadron- nucleus interactions. 

On the theoretical side, of great importance is the pQCD treatment of non-linear 
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interaction effects, in particular, concerning particle production in the transition 
region close to the parton saturation regime. Here, one may expect further progress 
to come from the "quasi-classical QCD" approach |29j . 

The author would like to acknowledge fruitful discussions with R. Engel. 
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